Density of states in d-wave superconductors of finite size 
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We consider the effect of the finite size in the a&-plane on the surface density of states (DoS) in 
clean d-wave superconductors. In the bulk, the DoS is gapless along the nodal directions, while 
the presence of a surface leads to formation of another type of the low-energy states, the midgap 
states with zero energy. We demonstrate that finiteness of the superconductor in one of dimensions 
provides the energy gap for all directions of quasiparticle motion except for = 45° (9 is the angle 
between the trajectory and the surface normal); then the angle-averaged DoS behaves linearly at 
small energies. This result is valid unless the crystal is 0°- or 45°-oriented (a 7^ 0° or 45°, where 
a is the angle between the a-axis and the surface normal). In the special case of a — 0°, the 
spectrum is gapped for all trajectories 9; the angle-averaged DoS is also gapped. In the special case 
of a — 45°, the spectrum is gapless for all trajectories 9; the angle- averaged DoS is then large at low 
energies. In all the cases, the angle-resolved DoS consists of energy bands that are formed similarly 
to the Kronig-Penney model. The analytical results are confirmed by a self-consistent numerical 
calculation. 

PACS numbers: 74.78.Bz, 74.78.Fk, 74.45.+C 



I. INTRODUCTION 

Among many possible types of unconventional 
superconductors^ the materials with the d-wave sym- 
metry of the pair potential are most widely dis- 
cussed. The d-wave symmetry is established in 
the quasi-two-dimensional high-temperature cuprate 
superconductors^ A characteristic property of the d- 
wave superconductivity is the gapless spectrum of quasi- 
particles. The pair potential is anisotropic, and the gap 
vanishes along the nodal directions. Another source of 
low-energy quasiparticles is the surface, which leads to 
forming the midgap states (MGS)j^ this phenomenon is 
due to the change of the sign of the pair potential along a 
trajectory upon reflection. Thus a d-wave superconduc- 
tor is a gapless superconductor with two types of low- 
energy quasiparticles. 

The d-wave superconductors can be employed in novel 
types of logic elements, qubits; 5 '- there is experimen- 
tal progress in this direction^S However, the low- 
energy quasiparticles introduce decoherence in d-wave 
qubits. 9,10 At the same time, the authors of Ref. men- 
tion a possibility to suppress the low-energy quasiparti- 
cles due to the finite size of the d-wave banks. 

In this paper, we study the influence of the finite size 
of a d-wave superconductor on the low-energy density of 
states (DoS) at the surface. There is only a limited num- 
ber of results related to this issue. The angle-averaged 
surface DoS was numerically studied by Nagato and Na- 
gai for 45°-oriented superconductors^ There is also a 
number of numerical results on the DoS in clean SN 
systems (where S is a conventional s-wave superconduc- 
tor and N is a normal metal), which can be relevant to 
the nodal directions of finite-size d-wave superconduc- 
tors due to similarity of the pair potential profile along 
quasiparticle trajectories in the two systems (see below 



for details). The works by van Gelder 12 and Gallagher 13 
are the most relevant in this respect. In Ref. ll4L She- 
lankov and Ozana suggested a method to treat multiple- 
interface superconducting systems and, as an application, 
numerically considered the DoS in a finite-size bilayer. 
Their results are relevant for the 45° trajectory in the d- 
wave system. Finally, we mention an analytical result of 
Ref. H5I where Fauchere et al. considered an SN system 
with repulsive interaction between the electrons in the N 
layer. In our language, their result refers to the splitting 
of the MGS. 

We demonstrate that the continuous normal-metallic 
spectrum along the nodal directions transforms to a set of 
energy bands; the lowest band is separated from zero by 
a gap. The midgap states (existing in the finite intervals 
of the trajectory angles 45° — a < \0\ < 45° + a - 
see Fig. n for notations) transform to two energy bands, 
situated symmetrically around zero, thus a gap appears. 
This gap depends on the direction of the quasiparticle 
trajectory, and vanishes for exactly = 45° direction (at 
any crystalline orientation a). After averaging over the 
directions, the DoS is zero at zero energy, and behaves 
linearly at small energies. Thus the averaged surface DoS 
remains gapless but is strongly suppressed compared to 
that in the bulk d-wave superconductor (i.e., at L — > 00). 

In Sec.m we present analytical results for the DoS. In 
Sec. IIIII we confirm and illustrate these results by self- 
consistent numerical calculations. Finally, we present our 
conclusions in Sec. IIVI 



II. ANALYTICAL RESULTS 



We consider a quasi-two-dimensional d^^^-wave su- 
perconductor of finite width, i.e., a strip in the aft-plane 
- see Fig. ^ A quasiparticle trajectory is sequentially 
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FIG. 1: d. 



-wave superconductor of finite width L in the 



afe-plane (quasi-two-dimensional strip). The orientation of 
the crystalline o-axis with respect to the surface normal is 
denoted by a, then the angular dependence of the pair po- 
tential is A(#) = Ao cos(2# — 2a). The pair potential along 
the trajectory described by the angle 6 changes periodically. 
We assume 0° ^ a 45° — this interval covers all physically 
different situations. 



reflected from one or the other surface of the strip, and 
the pair potential felt by the quasiparticle changes peri- 
odically. The profile of the pair potential along the tra- 
jectory is schematically depicted in Fig. [2] This profile is 
not self-consistent, while the self-consistent pair potential 
is suppressed near the surfaces (see Sec. lIIII below'). How- 
ever, the width of the regions where this happens has the 
characteristic scale of the coherence length £ = Vf/2ttT c 
(vf is the Fermi velocity, T c is the superconducting crit- 
ical temperature^). We assume L ^> £, then the regions 
where A is suppressed are relatively narrow, hence the 
piecewise constant A, depicted in Fig. [21 is a good ap- 
proximation. The results of self-consistent numerical cal- 
culations will be discussed below in Sec. IIIII they agree 
with the analytical results of the present section. 

The real-energy Eilenberger equations along the tra- 
jectory have the form (we choose A to be real)i£ 



-2iE + v F -^)j -2A/y = 0. 



-2iE-v F — j f-2Ag = 0, 
v F JLg + A(f-f)=0, 



(1) 



with the normalization condition 
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// = !■ 



(2) 



Here g is the normal Green function, while / and / are 
the anomalous Green functions that describe the super- 
conductivity. 



FIG. 2: The pair potential profile along a trajectory in the 
d-wave superconductor of finite width. The signs are chosen 
in such a way that the midgap states exist at Ai,A2 > 0. 
For definiteness, we choose Ai ^ A2; this does not impose 
any restrictions on the results. The pair potential changes on 
the surfaces of the superconductor. With respect to Fig. 
the introduced notations are defined as Ai = A(8), A2 = 
-A(-O), and d = L/ cos0. 



Since the Green functions along the trajectory are con- 
tinuous upon specular reflection at the surfaces of the 
d-wave superconductor, they obey the same condition 
in the effective one-dimensional problem described by 
Fig- [21 5; /) an d / must be continuous at the edges 
of the intervals of constant A. 

The DoS (normalized by its normal- metal value) is de- 
termined by the real part of the normal Green function: 



v = Keg. 

The DoS is symmetric, v{E) = v(—E). 
In the bulk, the Green functions are 



-iE 



<J 



VA 2 - E 2 ' 
which yields the standard DoS 

v = Re ■ 



/ = / 



VA 2 - £ 2 



\E\ 



(3) 



(4) 



(5) 



VE 2 - A 2 

The general solution of the Eilenberger equations (JTJ 




A 



G 



— I iE + VA 2 - E 2 
2 1 iE-VA' 2 -E 2 J 

A 



.kx 



with 



lE-^fW~^W_ 
JE + V& 2 -E 2 / 

, 2VA 2 - E 2 



-kx 



vf 



(6) 



(7) 
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The normalization condition J5J yields: 



When the pair potential is spatially symmetric, A(x) = 



A(— x), the Green functions obey the relations g(x) — 
g(—x), f(x) — f(—x), hence B = C. In our case, the 
pair potential is symmetric with respect to the center of 
each interval, thus we can write 



at — d < x < 0, and 



-iE 
Ai 



-iE 
-A 2 



Si 



Bo 



iE cosh(ki(x - 



Ai cosh(fci(a; 
d/2)) + VAf 



■d/2)) 

E 2 sinh(fci(x 



•d/2)) 



•d/2)) 



— A 2 cosh(fc 2 (ir 
i£cosh(fc 2 (;E - d/2)) + y/Aj - E 2 sinh(fc 2 (x - d/2)) 



(9) 



(10) 



at < x < d, while the solution on all the other intervals is obtained due to 2<i-periodicity of the Green functions. 

The boundary conditions (continuity of the Green functions) at x = and x = d yield four equations, only three of 
which are independent. Solving this system of equations and using the normalization condition JHJl, we determine all 
the four coefficients A\, A2, B\, and -B 2 . Then we find the Green functions and the DoS at the interface between the 
two intervals (this corresponds to the surface DoS in the <i-wave superconductor): 



v [x = 0) = Re ■ 



^A 2 - E 2 tanh 



E 2 tanh (Y) 



cosh(k\d/2) coshf/c2^/2) 



(ii) 



(E 2 + AiA 2 ) tanh (^) tanh ( - y/ A\ - E 2 ^/ A\ - E 2 



A similar formula was obtained by Gallagher at the 
same time, our results are different since his formula 
refers to the center of the strip in the d-wave problem, 
while we study the surface DoS. 

The zero-energy DoS can be found immediately. If 
Ai ^ A 2 , then Eq. mi yields v{E = 0) = 0. If Ai = A 2 
(= A), then we obtain 

v(E = 0) = cosh (^-^j ■ (12) 

In the above calculations, we did not assume Ad/vp 3> 
1 , however it is necessary to assume this in order to make 
the results physically sound, because in a d-wave su- 
perconductor, the piecewise constant pair potential (see 
Fig. |2) is a good approximation only under this condi- 
tion. 

Below we analyze Eq. at low energies, E <C A, 

in the following relevant cases: a) A 2 = (nodal direc- 
tions), b) Ai 7^ A 2 , and c) Ai = A 2 (the latter two 
cases correspond to the presence of the MGS in the infi- 
nite system). 

A. Effect of the finite size on the nodal 
quasiparticles 

Let us consider a nodal direction, which corresponds 
to A 2 = 0. For brevity, we shall denote Ai by A. 



In the bulk, the DoS along a nodal direction is normal- 
metallic, 

Vd^oo = 1- (13) 



The finite-size problem was considered previously (al- 
though in a different context) by van Gelder— and 
Gallagher. 13 While van Gelder considered an academic 
one-dimensional model with periodic step function A(x), 
in the work by Gallagher this situation emerged along 
a quasiparticle trajectory in a superconductor-normal- 
metal sandwich (with conventional s-wave superconduc- 
tor). They numerically demonstrated that in this su- 
perconducting version of the Kronig Penney model, 19 
the quasiparticle spectrum consists of energy bands with 
square-root singularities at the band edges. Below we 
present analytical results for this problem. 

Taking into account that 



£<A, (15) 
we obtain from Eq. I|ll|) : 



The bands are situated around the energies at which _Esin (Ed/vp) — A cos (Ed/vp) — 0. Since we consider E <C A, 
then cos(Ed/vp) must be very small to fulfill this equation, hence its argument is very close to ir/2 (for the lowest 
band) . Thus the center of the lowest band is 



7T v F f v F \ irv F 
~ 2~d V ~ Ad) ~ 2~d 



(17) 



In the vicinity of Eq , the DoS can be written as 



v — Re 



vpjd 



E-E Q + 2?fexp (- 



Ad 



E - E - 2^f exp 



(-*) 



(18) 



hence the width of the band is 
SE 



a Vf 

4 — exp 
d 



Ad 

Vp 



(19) 



and the DoS has integrable square-root singularities at 
the edges of the band. 

The minimal value of the DoS in the band is achieved 
at E — Eq, and equals 



1 Ad 

t'min = - exp 

2 \VF 



(20) 



The physical mechanisms behind the above results are 
quite transparent. Instead of the normal- metallic situa- 
tion that takes place for the nodal directions in the bulk, 
in the finite system we obtain the profile of the pair po- 
tential corresponding to the SN superlattice (Fig. [2] with 
A2 = 0). Then the energy spectrum in each normal layer 
consists of the Andreev levels^ which are smeared into 
the bands due to periodicity of the system. The energy 
of the lowest Andreev level corresponds to the center of 
the band, see Eq. (|17fl . while smearing is due to tunnel- 
ing across the barrier of height A and width d, and thus 
contains the tunneling exponential, see Eq. (|19|) . 



B. Effect of the finite size on the midgap states 

In the infinite system (L, d — > 00), the midgap states 
arise if the pair potential changes its sing upon reflection 



from the surfaced According to our definitions (Figs. 
and[5J, this happens at Ai, A 2 > 0. The MGS are local- 
ized near the surfaces and have exactly zero energy; the 
corresponding DoS is 



Vd- 



= 2tt- 



AiA 2 



-5(E). 



(21) 



Ai + A 2 

Below we consider the effect of finite d on this result; the 
results for the cases of differing and coinciding A x and 
A 2 will be qualitatively different. 



The case Ai > A 2 



Let us consider the general case, when Ai and A 2 are 
nonzero and can be different (still, the sings are such that 
the MGS exist in the infinite system). Employing 



Aid A 2 d 



» 1, 



vf vf 
E<A 1 ,A 2 , 



and expanding Eq. (|ll|l . we obtain: 



(22) 
(23) 
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This yields two bands, situated symmetrically around 
zero energy. The edges of the positive-energy band are 



E, 



2AiA 2 



6. min — 



2A X A 2 
Ai + A 2 



exp 



exp 



A 2 d 

VF 



A 2 d 

VF 



exp 



exp 



Aid 

VF 



Aid 

VF 



(25) 



(26) 

The minimal value of the DoS in the band is achieved at 



E 



2Aj A; 

Ai + A 



L V e 2 + e i> an( i equals 



1 1 



1 1 



= n1 1 + "2 ~ ^exp 



Aid 



(27) 



At the edges of the bands, the DoS has integrable square- 
root singularities. 

The position of the center of the band in the limit 
Ai ^> A 2 was calculated in Ref. fl5l (although in a differ- 
ent context), while only discrete energy levels were dis- 
cussed and the width of the band was not studied. 

Physically, the obtained results can be explained as 
follows. In the infinite system, d — > 00, we have the 
zero-energy levels localized near the interfaces between 
Ai and — A 2 . The first effect of the finite d is to split the 
levels at the two neighboring interfaces due to tunneling 
across the A2 barrier (the lowest barrier). This splitting 
is symmetric with respect to zero energy. Physically, it is 
similar to the standard quantum-mechanical problem of 
the level splitting in the double- well potential^ The sec- 
ond effect of the finite d is to smear each of the split levels 
due to periodicity of the system; the smearing is due to 
tunneling across the Ai barriers. This is similar to the 
standard Kronig-Penney models As a result, the center 
of the band (|25|) - (|26[1 is determined by the tunneling ex- 
ponential containing A2, while the width of the band is 
determined by the tunneling exponential containing Ai. 
Since A2 < Ai, the splitting of the zero-energy level is 
larger than its smearing, hence a gap in the spectrum 
arises. 



2. The case Ai = A 2 
At Ai = A 2 (= A), the DoS reduces to 



A / Ad 
- £ b = 2Aexp 



(28) 



Thus the MGS is smeared into the band of width 2Ej, 
around zero. The minimal value of the DoS is achieved 
at E = 0, and equals 



1 Ad 

"min = - exp 

2 \VF 



(29) 



while at the edges of the band (at E — ±-E&), the DoS 
has integrable square-root singularities. 

This result means that the two bands (positive and 
negative) that existed at Ai > A2, touch each other at 
E = and merge into a single band, while the singular- 
ities at E = transform into the minimum (|29|l . The 
equivalent result about the band centered at E = was 
numerically obtained in Ref. 14 (although in a different 
context). 

The physical explanation of these results is the same as 
in the previous case, Ai > A2. The only difference is that 
in the case of equal barriers, Ai = A2, the splitting of 
the zero-energy level is exactly the same as its smearing, 
hence no gap in the spectrum appears. 



C. Angle-averaged DoS 

Let us consider the behavior of the angle-averaged sur- 
face DoS at E — > 0. The only contribution to the DoS 
arises from the vicinity of the angles at which Ai = A2 — 
these are always the 8 — ±45° angles [at any orientation 
a, since A(0) = A cos(20 - 2a)]. 

Introducing ?9 via = tt/A — we expand the pair 
potential in the vicinity of 6 = 45°: 

Ai w A + A'tf , A 2 w A - A'0, (30) 

where 

A = A sin 2a, A' = 2A cos 2a. (31) 
Then we expand 



e-i — ei w exp 



we assume that 



VF 



Ad\ 2A'M 

V F ) V F 



< 1 



(32) 



(33) 



(since A' ~ A, this condition implies that d -C Up/ Ad 
1). The angle-resolved DoS 124|) then takes the form 



6 



= A Re 



E 




-{E 2 - 



2AA'tfd 



: exp 



(-£) 



E 2 - 



2Aexp 



Ad 



(34) 



The contribution to the angle-averaged DoS comes from 
the interval — < i9 < t9j,, where 



2AAVi _Ad 

Hi? ~ V 



(35) 



Condition 13311 for all $ up to fib is equivalent to 

E<t:E b , (36) 

where £k = 2A exp (— Ad/vp) is the upper edge of the 
band. 

The angle-resolved DoS (|34|l can be written as 



1 



4AA 



Re- 



E 



^exp(-^) VWE) 



.02 



(37) 



The angle-averaged DoS is 



E 



2AA 



^exp(-^ 



(38) 



The DoS is zero at zero energy, and behaves linearly at 
small E. Thus the averaged surface DoS remains gapless 
but is strongly suppressed compared to that in the bulk 
d-wave superconductor (i.e., at L — > oo). 

The gapless structure of the angle-averaged DoS 
is due to integrating in the vicinity of 9 = ±45°, where 
the gap in the angle-resolved DoS vanishes. At the same 
time, the DoS probed by transport methods^, is given by 



ftr 



v{9)D{9)d9, 



(39) 



and differs from Eq. lj3*5|) by the weighting factor D(9), 
the angle-dependent transparency of the tunneling in- 
terface. When the tunneling interface has finite thick- 
ness, the function D{9) is exponentially suppressed at 
not too small 9, then the contribution of the 9 = ±45° 
trajectories can be significantly suppressed [since v tI ~ 
D(45°)^av at low energies]. In this case, we can ex- 
pect that the transport DoS vtr will be gapped despite 
the angular averaging. The directional selectivity of 
tunneling is most pronounced in the scanning tunnel- 
ing spectroscopy experiments, where the effective tun- 
neling cone around the surface normal can be as narrow 
as 59 ~ 20° ^ 

The result l|38|l for the angle-averaged DoS does not 
refer to the cases a = 0° and a = 45° (and close orienta- 
tions); these cases are special. 



At a = 0°, the MGS do not appear, and the low-energy 
DoS is entirely due to the nodal directions. Then accord- 
ing to Sec. Ill Al the spectrum along a nodal direction 
acquires a gap due to the finite size. In this situation, 
the angular averaging preserves the gap in the spectrum, 
approximately given by Eq. i|17|) ■ 

At a = 45°, the condition Ai = A2 (which implies the 
gapless spectrum) is satisfied not only at 9 — ±45° but 
at any 9 (all trajectories). Then angular averaging does 
not introduce new qualitative features (compared to the 
angle-resolved result), and the DoS at small energies is 
large according to the results of Sec. Ill B 21 Comparing 
with the bulk case, we can say that due to the finite 
size, the zero-energy peak in the DoS is smeared but not 
split. This statement agrees with the results of Nagato 
and Nagaiil who studied this case by a self-consistent 
numerical method. 

The above results refer to the surface DoS. However, 
the averaged DoS is linear at low energies also inside 
the strip, although the slope is different from Eq. 138|l . 
because the MGS contributing to this result decay into 
the bulk of the sample. For example, in the middle of 
the strip the angle-averaged DoS differs from Eq. 138|) 
by an additional factor 2exp(— Ad/vF )■ This result is 
again valid at E <C E^; this interval shrinks in the limit 
of large d, where the DoS is mainly determined by the 
standard nodal contribution at larger E. 



III. SELF-CONSISTENT NUMERICAL 
RESULTS 

In order to take into account the spatial inhomogeneity 
of the pair potential, we solve the problem numerically 
by a self-consistent numerical method similar to the one 
used in Ref. |3 We rewrite the Eilenberger equations 
employing the Riccati parametrization^ which is well 
suited for numerical integration. Starting from the spa- 
tially homogeneous pair potential, we calculate the Green 
functions over the whole sample, and then use the self- 
consistency equation to refine the initial approximation 
for the pair potential. The procedure is repeated itera- 
tively until the necessary accuracy is reached. 

An example of the pair potential along a trajec- 
tory, calculated self-consistently, is shown in Fig. [3J In 
anisotropic superconductors, surfaces lead to pair break- 
ing and suppression of the pair potential. The region 
where the suppression takes place, has the characteristic 
scale of £, hence at L ^> £ the piecewise constant A (see 
Fig. [2| is a good approximation. 



7 




0.4- 



0.2 H , , , , . 1 . 1 

-40 -20 20 40 

FIG. 3: Self-consistent profile of the pair potential along the 
8 = 40° trajectory at the crystalline orientation a = 20°. 
The width of the rf-wave strip is L = 20 £, the length of the 
trajectory between the successive reflections from the surfaces 
is d = Lj cos 40°. 
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FIG. 4: The low-energy angle-resolved DoS for several trajec- 
tories 9 at the width of the strip L = 10 £ and the crystalline 
orientation a = 20°. 

Figures 01 and [5] demonstrate the angle-resolved DoS 
for two values of the strip width: L = 10 £ and L = 
20 £. These numerical results agree with the analytical 
theory developed in Sec. [H] Indeed, we observe that the 
spectrum consists of the energy bands, both along the 
nodal direction (9 = 25°)S& and the directions for which 
the MGS exist in the limit L -+ oo (9 = 30°, 35°, 40°, and 
45°). The DoS is minimal in the centers of the bands and 
demonstrates the square-root singularities on their edges. 

The energy gap for the lowest band along the nodal 
direction (9 = 25°) is larger than for the MGS direc- 
tions; this agrees with our analytical results according to 
which the gap for the nodal direction does not contain 
an exponentially small factor [see Eq. (|17|l and note that 
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FIG. 5: The same as in Fig. H but for L = 20 £. 



the energy gap is of the same order as the center of the 
band since the width of the band is much smaller] . The 
lowest energy band for the nodal direction can be seen 
in Fig. while in Fig. [5] it is beyond the demonstrated 
energy range. 

The energy bands for the MGS directions are gapped if 
9 45°, which corresponds to the case Ai ^ A2 studied 
in Sec. Ill B II The 9 = 45° direction corresponds to the 
case Ax = A 2 studied in Sec. Ill B 21 the energy band is 
gapless in this case, its center is E = where the DoS 
reaches a minimum. 

As L increases, the energy gaps become smaller and 
the bands become narrower. 

The numerical results allow us to directly check the 
accuracy of our analytical estimates. For example, we 
consider the 9 = 40° direction at L — 20 £. Then the 
edges of the energy band are given by Eqs. (|25|l and 
P| which yield the center of the band E /T c = 0.015 
and the width of the band SE/T C = 0.0029, while the 
numerical results are E /T c = 0.017 and SE/T C = 0.0037 
(see Fig- - Thus we conclude that our analytical results 
provide reasonable quantitative accuracy while capturing 
all qualitative features of the DoS. 



IV. CONCLUSIONS 

We have considered the effect of the finite size in the 
a6-plane on the density of states (DoS) in d-wave super- 
conductors. We assumed finiteness in one direction, thus 
we deal with quasi-two-dimensional strip instead of quasi- 
two-dimensional plane. We assumed arbitrary crystalline 
orientation, described by the angle a between the a-axis 
and the normal to the surfaces. 

The problem is solved analytically neglecting the sup- 
pression of the pair potential near the surfaces; the results 
are confirmed by a self-consistent numerical calculation. 
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In the relevant limits, our results agree with previous 
studies. 

In the bulk, the DoS is gapless along the nodal direc- 
tions, while the presence of a surface leads to formation of 
another type of the low-energy states, the midgap states 
with zero energy. Due to the finite size of the super- 
conductor, the spectrum of nodal quasiparticles acquires 
an energy gap. The midgap states acquire the angle- 
dependent gap that vanishes at 9 — 45° (9 is the angle 
between the trajectory and the normal to the surface); 
this result is valid unless the crystal is 0°- or 45°-oriented 
(a 7^ 0° or 45°). In the special case of a — 0°, the midgap 
states are absent, and the spectrum is gapped for all tra- 
jectories (including the nodal directions). On the oppo- 
site, in the case of a = 45°, the spectrum is gapless for 
all trajectories 9. In all the cases, the angle-resolved DoS 
consists of energy bands that are formed similarly to the 
Kronig-Penney model. 

In the special case of a = 0°, the angle-averaged DoS 
has a gap. In the special case of a — 45°, the angle- 
averaged DoS is finite at low energies. 

At a ^ 0° or 45°, the angle-averaged surface DoS is 
strongly suppressed due to finite size, while remains gap- 
less and behaves linearly at small energies. The low- 
energy contribution comes from the trajectories with 
9 ps 45°, hence we can expect that the energy gap sur- 
vives upon angular averaging if one measures the trans- 
port DoS in the case when the 45°-angle contribution is 
suppressed by the transparency of the tunneling barrier. 

In our model, we did not take into account such im- 
perfections in the system as bulk impurities2£i22i22i22i2i 
and interface roughness iSSM We expect them to smear 
the features discussed in the present paper. In particu- 
lar, the energy bands (in the angle-resolved DoS) sepa- 



rated from zero by gaps, exist only when imperfections 
are weak, otherwise the bands are smeared and finite DoS 
appears at zero energy. A rough estimate demonstrates 
that smearing of the bands due to the bulk impurities ex- 
ceeds the width of the bands if the electronic mean free 
path I is smaller than the width of the strip L. However, 
the effects discussed above can survive in the case of rare 
impurities, when I 3> L and long intervals along the strip 
(between two neighboring impurities) can be considered 
clean. 

Our model is not material-specific and can be in prin- 
ciple applied to any rf-wave superconductor. However, 
the only class of materials where the d-wave symmetry 
is established at present, is the high-T c oxides, where the 
coherence length £ is very small. Then the finite-size ef- 
fects discussed above become pronounced at quite small 
size L. 



Acknowledgments 

We are grateful to M. H. S. Amin, Yu. S. Barash, 
A. M. Bobkov, I. V. Bobkova, M. V. Feigel'man, 
M. Yu. Kupriyanov, G. Leibovitch, and A. M. Zagoskin 
for helpful discussions. The research was supported by 
the D-Wave Systems Inc. and the ESF PiShift program. 
Ya.V.F. was also supported by the RFBR grants Nos. 03- 
02-16677 and 04-02-16348, the Russian Science Support 
Foundation, the Russian Ministry of Industry, Science 
and Technology, the program "Quantum Macrophysics" 
of the Russian Academy of Sciences, CRDF, and the Rus- 
sian Ministry of Education. 



* Electronic address: fominov@landau.ac.ru 

' Electronic address: a.golubov@tnw.utwent e"nl| 

1 V. P. Mineev and K. V. Samokhin, Introduction to Un- 
conventional Superconductivity (Gordon & Breach, Ams- 
terdam, 1999). 

2 D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 

3 C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

4 C.-R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 

5 L. B. Ioffe, V. B. Geshkenbein, M. V. Feigel'man, 
A. L. Fauchere, and G. Blatter, Nature (London) 398, 679 
(1999). 

6 A.M. Zagoskin, |cond""m at/9903170 (1999); A. Blais and 
A.M. Zagoskin, Phys. Rev. A 61, 042308 (2000). 

7 E. Il'ichev, M. Grajcar, R. Hlubina, R. P. J. IJs- 
selsteijn, H. E. Hoenig, H.-G. Meyer, A. Golubov, 
M. H. S. Amin, A. M. Zagoskin, A. N. Omelyanchouk, 
and M. Yu. Kupriyanov, Phys. Rev. Lett. 86, 5369 (2001). 

8 A. Ya. Tzalenchuk, T. Lindstrom, S. A. Charlebois, 
E. A. Stepantsov, Z. Ivanov, A. M. Zagoskin, Phys. Rev. B 
68, 100501(R) (2003). 

9 Ya. V. Fominov, A. A. Golubov, and M. Yu. Kupriyanov, 



Pis'ma Zh. Eksp. Teor. Fiz. 77, 691 (2003) [JETP Lett. 
77, 587 (2003)]. 

M. H. S. Amin and A. Yu. Smirnov, Phys. Rev. Lett. 92, 
017001 (2004). 

1 Y. Nagato and K. Nagai, Phys. Rev. B 51, 16254 (1995). 

2 A. P. van Gelder, Phys. Rev. 181, 787 (1969). 

3 W. J. Gallagher, Phys. Rev. B 22, 1233 (1980). 

4 A. Shelankov and M. Ozana, Phys. Rev. B 61, 7077 (2000). 

5 A. L. Fauchere, W. Belzig, and G. Blatter, Phys. Rev. Lett. 
82, 3336 (1999). 

6 The zero-temperature amplitude of the pair potential in 
the bulk is related to the critical temperature as Ao = 
2.14T C , where the amplitude Ao is defined so that A(8) = 
A o cos(20 - 2a). 

7 G. Eilenberger, Z. Phys. 214, 195 (1968). 

8 It is interesting to note that Gallagher solved the Gor'kov 
equations with subsequent averaging over spatial oscilla- 
tions on the scale of I/Zcf (&f is the Fermi wavevector), 
while we start from the quasiclassical Eilenberger equa- 
tions. Thus the correspondence between our results im- 
plies that the quasiclassical approximation in the problem 
is valid provided IzfL ^> 1. 



9 



R. de L. Kronig and W. Penney, Proc. R. Soc. London, 
Ser. A 130, 499 (1931). See also S. Fliigge, Practical Quan- 
tum Mechanics I ( Springer- Verlag, Berlin, 1971). 
A. F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964); 49, 
655 (1965) [Sov. Phys. JETP 19, 1228 (1964); 22, 455 
(1966)]; I. O. Kulik, Zh. Eksp. Teor. Fiz. 57, 1745 (1969) 
[Sov. Phys. JETP 30, 944 (1970)]. 

L. D. Landau and E. M. Lifshitz, Quantum Mechan- 
ics: Non-Relativistic Theory, Volume 3 (Butterworth- 
Heinemann, Oxford, 1981). 

Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, 3451 
(1995). 

A. Sharoni, G. Leibovitch, A. Kohen, R. Beck, 
G. Deutscher, G. Koren, and O. Millo, Europhys. Lett. 
62, 883 (2003). 

M. H. S. Amin, A. N. Omelyanchouk, S. N. Rashkeev, 
M. Corny, and A. M. Zagoskin, Physica B 318, 162 (2002). 
N. Schopohl and K. Maki, Phys. Rev. B 52, 490 (1995); 
N. Schopohl, cond-mat/9804064 (1998). 
At a = 20°, one of the nodal directions in the bulk is 9 — 
—25°. In the strip geometry, we call the 6 = 25° direction 



nodal because the corresponding trajectory includes the 
nodal 6 — —25° segments (see Fig. 0. 
P. J. Hirschfeld, P. Woffle, and D. Einzel, Phys. Rev. B 
37, 83 (1988). 

Yu. S. Barash, A. A. Svidzinskii, and V. P. Mineev, Pis'ma 
Zh. Eksp. Teor. Fiz. 65, 606 (1997) [JETP Lett. 65, 638 
(1997)]. 

A. Poenicke, Yu. S. Barash, C. Bruder, and V. Istyukov, 
Phys. Rev. B 59, 7102 (1999). 

M. L. Kulic and O. V. Dolgov, Phys. Rev. B 60, 13062 
(1999). 

Y. Asa no, Y. Tanaka, and S. Kashiwaya, 
cond-mat/0302287 (2003). 

M. Matsumoto and H. Shiba, J. Phys. Soc. Jpn. 64, 1703 
(1995). 

Yu. S. Barash, A. A. Svidzinsky, and H. Burkhardt, Phys. 
Rev. B 55, 15282 (1997). 

A. A. Golubov and M. Yu. Kupriyanov, Pis'ma Zh. Eksp. 
Teor. Fiz. 69, 242 (1999) [JETP Lett. 69, 262 (1999)]. 



